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Abstract. We give a CLR type bound on the number of bound states of Schro- 
dinger operators with matrix-valued potentials using the functional integral method 
of Lieb. This significantly improves the constant in this inequality obtained earlier 
by Hundertmark. 

o 

! 1. Introduction 

We consider the Schrddinger operator —A — V(x) on M. d , but with the difference 
from the usual case that V is a Hermitian matrix-valued potential. In other words, 
the Hilbert space is not L 2 (IR d ) but L 2 (IR d ; C N ). The values of functions in this space, 
i/j(x), are N— dimensional vectors. (What we say here easily generalizes to 'operator- 
valued' potentials, i.e., is replaced by a Hilbert space such as L 2 (M m ), but we stay 
t^. ; with matrices in order to avoid technicalities.) The Cwikel-Lieb-Rozenblum (CLR) 

bound for d > 3 in the scalar case N — 1 states that #(— A — V), the number of 
negative eigenvalues of —A — V, can be estimated by 

#(- A -V)< L 04 [ V + (x) d ' 2 dx . (1.1) 



o 



(Here and below v± := (\v\ ± v)/2 denotes the positive and negative part of v.) We 
remind the reader that the 'semi-classical' approximation to #(— A — V) is given in 



5_i ■ the scalar case by the phase space volume 



(2n)- d [[ dpdx = L c l d [ V+{x) d/2 dx 



where 



Lt 4 = W d I d P= (2V/ 2 r(d/2 + l))" 1 

J{p&R d : p 2 <l} 
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The bound (11.11) was obtained by completely independent methods in [CllLlEj- Later, 
different proofs were given in \Co\ ILY] . The best constant, which is close to optimal 
for d = 3, was obtained in [LJ using the Feynman-Kac formula and Jensen's inequality. 

Our goal here is to extend inequality ( 11.11) to the matrix case (with a possibly 
different constant Lo,d) • The motivation for this extension was the work of Laptev and 
Weidl |LWlj who realized that the extension allowed one to conclude that good/sharp 
constants obtained in low dimensions would automatically give good/sharp constants 
in higher dimensions. The fact that the inequality (11.11) is valid in the matrix case was 
proved by Hundertmark [H] , confirming a conjecture in [LW2] . He follows Cwikel's 
method and obtains a constant which is far from optimal. Hundertmark points out 
that 'it would be nice to extend Lieb's [. . . ] proof of the CLR-bound to operator-valued 
potentials'. This is the content of this letter. 

Theorem 1.1. Let d > 3 and assume that V is a function on M. d taking values in the 
Hermitian N x N matrices. Then 

#(-A -V)< Ro, d L« d [ Tr cN [V + (x) d ' 2 ] dx (1.2) 

where R 0td < 10.332 and V+ := (\V\ + V)/2. 

The constant 10.332 will be obtained for d = 3 and, by the Laptev-Weidl method 
(as used by Hundertmark [H]) it is valid uniformly for all d > 3. We emphasize that 
our bound on i? ,d is slightly worse than the constant 6.87 in [L] for the scalar case 
N = 1. Still, it improves that of [H] by almost one order of magnitude. For d = 3 our 
bound on i? ,3 is a t most a factor 2.24 bigger than the optimal constant in (11.21) . since 
it is known that R , 3 > 8/V3 « 4.619 [LTj . 

It is well known that by a simple integration the bound ( 11. 2\\ yields the Lieb-Thirring 
inequalities 

Tr La(lW (-A - V)l < Ry, d L« [ Tr cN [V + { X y +d / 2 ] dx (1.3) 

for all 7 > 0, d > 3 with R l4 < R 04 < 10.332 and 

L* d = (2n)- d I (l-p 2 )ldp. (1.4) 

Indeed, R l4 is a monotone non- increasing function of 7 [AL] . Even in the scalar 
case N — 1, this yields the best known constants in this inequality for the parameter 
range < 7 < 1/2. For comparison we recall that the best known bounds for larger 
values of 7 are R l4 < 2n/V3 w 3.628 if 7 > 1/2 and R l4 < tt/VZ « 1.814 if 7 > 1 
jHLWl [DLL] . For 7 > 3/2 one has Ry >d = 1, which is sharp [LWl] . We refer to the 
surveys [HI ILW2j for more about inequalities (11. 3ft . 

Apart from yielding very accurate constants we believe that there is a mathematical 
interest in extending the path-integral method in [L] to the operator-valued situation. 
In contrast to the method of [C] used in [Hj. which is rather rigidly based on mapping 
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properties of the Fourier transform, the method of [L] used here works in much wider 
generality, e.g. on Riemannian manifolds. The only input needed is an upper bound 
on the heat kernel of the (scalar) unperturbed operator. For example, the Hardy-Lieb- 
Thirring bounds in |FLS] extend to the matrix-valued situation. 

As already pointed out, we proceed similarly to [L]. Therefore we will be brief at 
some points and ignore some technicalities. There is an important new ingredient in 
our proof, however. Since matrices Wi, . . . , W n do not commute, in general, we need to 
work with the "time ordering" of a function f(Ylj Wj) °f their sum. In Proposition 13. II 
we shall prove a modification of Jensen's inequality valid in this setting for a certain 
class of convex functions /. 



2. A TRACE FORMULA 

Given self-adjoint iV x A-matrices W\, . . . , W n and a function / on R, the usual ma- 
trix f(J2j Wj) is defined by the spectral projections of J2j Wj. Instead, we introduce 
the "time-ordering" of the matrix f(^2j Wj) as follows. We write Wj in its spectral 
representation 

N 

k=l 

where Wj? are the eigenvalues and Pff the corresponding orthogonal projections, and 
define 

N / n \ 

Tf( Wl ,...,W n ):= £ f ^•••< n) " (2- 1 ) 

fci,...,fe„=i \i=i / 

Intuitively, this means that when calculating f(%2j Wj), one puts all the Wis left of 
the W 2 : s, the M^'s left of the W^s, and so on, without worrying about commutators. 
It is instructive to look at some examples. 

Example 2.1. If /(//) = fi k , k G N, then the definition immediately implies 

Tf{Wt,...,W n )= ' • . WF'-'W*. 

n+---+jn=k 

Example 2.2. If /(//) = e a ^, a6l, then again by the definition (12. ip 

Tf(W 1 ,...,W n )=e aWl ---e aW *. 
Similarly, one shows that if f(fi) = fie atl , a6l, then 
Tf(W x ,...,W n ) 

= W x e aWx e aW * ■ ■ ■ e aWn + e aWl W 2 e aW2 ■ ■ ■ e aWn + ... + e aWl e aW2 ■ ■ ■ W n e aWn . 



We have introduced the notion of time-ordering in order to generalize the trace 
formula in [L], which is the starting point of the analysis leading to (11. ip . 
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Proposition 2.3. Let f be a non-negative, lower semi- continuous function f with 
f(0) = 0, and let 



F(X) :-. 



poo 

/ /(/i)e^ /A /i -1 ^, A>0. 
Jo 



(2.2) 



Then for any sufficiently regular and decaying functions V on IR d ; d > 3, taking values 
in the non-negative N x N- matrices, one has 



Tr 



L 2 (M. d ;C N ) 

°° dt 



F(v 1/2 (-/\y l v l/2 ) 







lim 

t n^oo 



dx~^ m doCfi 



(2.3) 



i=i 



Tf -V( Xl ) 
n 



n 



-V{x r , 



with the convention that 



In the limit n — > oo the multiple integral on the right side of (12. 3p converges to 
a Wiener integral (the Feynman-Kac integral); in fact, the right side of ( 12.31) is the 
Trotter product approximation to this integral (TJ IRS} IS2] . 

Proof. By an approximation argument [S2, Thm. 8.2] it suffices to prove this formula 
for 

F(A) = A/(l + aA), f{ f i) = f ie- a », 

where a > is a constant. Using the resolvent identity and Trotter's product formula, 
one easily verifies that in this case 

F(V 1 /2(_A)- 1 V 1 / 2 ) = V^ 2 (-A + aVy l V 1 ' 2 



/■oo 

/ V l/2 exp(-t(-A + aV))V 1/2 dt 
Jo 

lim T n (t) dt . 



Here, 

T n (t) := V 1/2 (exp(tA/n)exp(-taV/n)) n V 1/2 . 

The latter is an integral operator and we evaluate its trace by integrating its kernel 
on the diagonal. Let k denote the heat kernel 

k(x,y,t) := (47Tt)- d / 2 exp(|x-y| 2 /(4t)). 

Then 

r ^ r L 2 (R d ;C N ) T n (t) 

t 



dxi ■ ■ ■ dx n TT k 



7? 



Tr f 



e -?^i)... e -?%) 
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Cyclical relabeling of the variables leads to 

(K d ;C w ) T n {t) 

J dxi ■ ■ ■ dx n Y\ k (xjjXj-x, Tr C iv [Tf(tV(x 1 )/n, . . .,tV(x n )/n)} 

(compare with Example 12.21) . The claimed formula ( 12.31) follows if one interchanges 
the trace with the t-integration and the n-limit. ■ 




3. Jensen's inequality and time ordering 

To apply (12.31) we need to estimate the trace of a time-ordered sum. Recall that 
Jensen's inequality says that Tr/Q^Wj) — n^ 1 J^Tr f(nWj) for / convex. The 
analog for the time-ordered case, and a certain class of /'s, is 

Proposition 3.1. Assume that 

CO „ 

/W = Vy+ / e-°*dn(a) (3.1) 

for some ao, a\ G a-,- > for j > 2 and a non-negative measure \i. Then for any 
non-negative N x N -matrices Wi, . . . , W n 

1 n 

Re Tr C iv [Tf(W 1} . . . , W n )\ < - £ Tr C iv f{nW 3 ) . 

H 3=1 

Note that the / in (13.11) is convex. We do not know whether the statement is true 
for an arbitrary convex function. If it were, the constant in Theorem 11.11 could be 
improved, as explained at the end of this letter. 

Proof. By linearity of the trace it suffices to consider the cases f\ji) = fi k , k G N, and 
/(//) = e QA \ In the former case, one has by Holder's inequality for traces (see, e.g., 
Thm. 2.8]) 

u 

Re Tr , x [Tf(W u . . . , W n )\ = £ rj — - Re Tr c ^ [W? ■ ■ ■ W£] 

jl+ ... +jn =k Jl -'" Jn - 

< £ ~~~\~~7~\ (Tr C iv Wi) jl/k ■ ■ ■ (Tr c jv W^) jn/k 



3l+-+3n=k 



f [ £ (Trc W f) 

d=i 



k\l/k 
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and the assertion follows from the convexity of /. In the latter case, one has similarly 
by Holder's inequality and the geometric-arithmetic mean inequality 

Re Tr c ^ [Tf(W 1} . . . , W n )] = Re Tr c ;v [e aWl ■ ■ ■ e aW "} 

< (Tr cN e anW ^) 1/n ■ ■ ■ (Tr cN e mW ») X,n 
1 " 

< - V Tr C iv e naW i , 
n / — ' 

3=1 

as claimed. ■ 

Corollary 3.2. Assume that f is a non-negative function of the form considered in 
Proposition \3. 1\ and let F be as in (12. 2p . Then for any sufficiently regular and decaying 
function V on IR d taking values in the non-negative N x N -matrices, one has 

TW»>f(V 1/2 (-A)-V^) < ^ (jf \V(*n d X . 

(3.2) 

Proof. Combining Proposition 13. II with Proposition 12.31 we obtain 

lV i2(Md;CW)J P(y 1 /2 ( _ A) -iyi/2 ) 

f°° dt f f / t\ 1 " 

< — lim / •■■/ Y\k I Xj,Xj-x, - ) -yjTVcJv f{W(x j )) dx x ■ ■ -dx n . 
Jo t ™y Rd Jteifi \ n) nj^ 

(Here we have used that the left side of (I2.3P is real, hence only the real part of Tr T f 
contributes to the integral.) The semi-group property implies 



- / ■••/ J7 k Xj-i, - ) ^ Trcw fjtVjxj)) dx T 



3=1 

= ~ / k *) Tl ':' A f( tV ( x j)) dx i = (A i d /2 / IV; ' V f( tV ( x )) dx ■ 

Denoting the eigenvalues of V(x) by Vi(x) < . . . < v N (x) one finds that 

/■°° dt Tr cj y f(tV(x)) ^ f°° dt f(t Vj (x)) A d/2 rdsf(s) 

Jo t t^ '^Jo t t^ j£ iKX) J 8**/*' 

thereby proving the assertion. ■ 

4. Proof of Theorem 11.11 

First we assume that d = 3. By the variational principle we can assume that V(x) 
is a non-negative matrix for all x, and by an approximation argument we can assume 
that V is smooth and rapidly decaying. For any increasing function F on (0, oo) the 
Birman-Schwinger principle implies that 

# (_ A _ V) < F (i)-i Tr i2(R3 . cW) F(y 1 /2 ( _ A) -iyi/2 ) (41) 
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We choose F = F a of the form (12.20 where a > is a parameter and / = f a is defined 
by 

u 2 a 2 r°° 

fju) = -J_ = p - a + = // - a + a 2 / e ~' (M+a) dt . 

fi + a p + a J 

Since this function is of the form considered in Proposition 13.11 we can apply Corol- 
lary 13.21 and get in view of (14. II) 

#(-A - V) < C a [ Tt c n [V(xf 2 ] dx 

where 




The result follows by choosing a = 1.13, which approximately minimizes C a . 

Now we assume that d > 4. We will use the Laptev- Weidl strategy to reduce this 
case to the case d = 3 as in [Hj. We note that by a straightforward approximation 
argument as in [LWTj the inequality for d = 3 holds also for N = oo, i.e., if V(x) 
assumes values in the compact self-adjoint operators on a separable Hilbert space. 
Introduce variables ) G M d where X\ G M 3 and x<i G M d 3 . We decompose 

the Laplacian correspondingly as —A = — Aj — A 2 and define, for fixed X\ G M 3 , 
W(a;i) := (—Ax — V(x\, •))_. If V is, say, smooth with compact support, then W^xi) 
is a compact operator in L 2 (IR d ~ 3 , C ) for every x\. The variational principle and the 
inequality for d = 3 imply that 

#(-A -V)< #(-Ai - W) < i2^ 3 / Tr i2(Rd - 3)CW) [^(xx) 3 / 2 ] dx x . 

By the result of Laptev and Weidl |LW1] , one has 

Tr i2(Rd -3 iCJV) [^(xx) 3 / 2 ] < L^ 2id _ 3 / Tr cN [V(x 1 , x 2 f 2 ] dx 2 

with the constant Ly 2d _ 3 from (jl.4p . Noting that -^0,3^3/2 d-3 = -^o,d we obtain the 
assertion of Theorem 11.11 

Remark 4.1. If the estimate in Proposition 13 . 1 1 held for all convex functions (not merely 
for those of the form (13. ip ). then we could choose / a (/x) = (// — a) + in the preceding 
proof, as in [L], and would get the same constant as in the scalar case. 
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